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MODIFIED MATRIX INVERSION - REAL NUMBERS 


DECUS Program Library Write-up 


DECUS NO* 8-550 


This is similar to "Matrix Inversion - Real Numbers" (DECUS NO. 8-72) by 
Professor A. E. Sapega. It has been modified to run under PS/8. Input is 
from the high speed paper tape reader and the output is routed to the DECwriter 
with the option of having a paper tape made of the inverse. 

















matrix INVERS I 


ON 


SIZE OF MATRix= 


VES = 0? L N O =i htR THPE 0UTPUT C-ESIRED? “- ^ -- 

™ of form, m cmim 



MATRIX INVERSION RESULTS 


original matrix values 


o. 200000002+01 


ROW l 

1y000000E+01 
ROW 2 

0 600000002+01 
ROW j: 


O- 79999999E+01 

the inverse matrix 


ROW 1 


0 50000000E+01 

0. 700000002+0^ 
IS 


0. 3000001002+01 
0-399999992+01 
0. S99999992+01 


"0. 80952380E-00 
ROW 2 

0- 104761902 + 01 
ROW 2 


~0. 14285 714E-00 
0- 71428572E-00 


0 333333332-00 
•0, 66666666E-00 


-0. 42857142E- 


00 


33333332E-00 


-0 95228089E-01 





























5 

103 


3 00 


3 ^ S' 


3B5 


3 06 


3 0? 


i o 

3 3 0 


WRITE <3, 3 3 
FORMfiT W' MfiTRIJ 
WRITE < 3 ,2 3 
F O ft M fi T < / 2 c, ] 7 £ 

J?EfiD. _<l, 2 3 n 
FORM fiT <12 3 
WRITE <3.. 62tO 
WRITE <3,6253 
REfif- <3,63:03 PLINTH 
WRITE <3,6253 ID, 
WRITE <3,6363 
PfioSF 
L R S 7 - N * N 
COMMON Fi 

DIMENSION FK 30003 
13-3 

I _2 = L fiS T ~ Nh 3 
f>0 3 00 J ~3,N 
00 3 01 1 = 11, 12, N 
REF;(3 <2,53 fi < I 3 
FORMfiT <E15. S3 
CONTINUE 
_ _ 11 = 11_+1 
12= 12+1 
CONTINUE 
WRITE <3,83 
FORMfiT </?/'‘ MfiTRlX 
WRITE <3,93 
FORMfiT <//'• ORIGINAL 
JOG = 3 
00 TO 200 
DO 3 50 J = l, n 
00 3 05 1=3, N 
fi<LfiST-t I 3 = 0. 0 
CONIINUE 
fi < L f; S 7 + J 3 = 3.0 
PVT=fl<J3 
J 2 = t fi ST +J 
00 3 06 KF'- Jj J2, N 
fi < K F’ 3 = ft < K E’ 3 / F’ V T 
CONTINUE 
DO 110 KRT*T,' n 
IF < KRT - J 3 10?, ai 0 , 107 
KR3=KRT 
KR2=KR1 + L RST 
K P R = J 

ft W C = fi < K R1 > 

DO 109^ KR - Kfti, KR2, N 

fi < K R 3 = fi < K R 3 - ft n C * fi < f p ft ) 

K. P ft = K ft R + N 

CON/iNUt 

CONiiNOE 

DO 111 j=i,LRST 


INVERSION 
OF M fi 7 RIX = 


3 


INVER 
M Fit ft 


SI 

IX 


ON 


RES 

flue 


Ut T 


S' 2/3 


3 
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11 a 

a by 

: 60 


.16 a 
200 


201 


>>< i - H <; a 4 h.) 

CUN I .1 NUt 
CUNIJ NUt 
WK'li 'fc <1, 16 0.) 

FORMAT <77 ? THE 

J OG = 2 

e jog=flWtTjog') 

GO TO 200 
CALL EXIT 
IL N E - <5 
NR -1 

NR' C - 1 _ 

WRI [ I a, 20 iy~NR~ 

FORMAT </'■ ROW ' , I, 
11-1 ~ — 
12 = L ft 5 T - N +1 


INVFRSE MATRIxT 



KONT = 1 “ .— — 

DO 210 JP=1,N 

JL~‘ 

DU 212' I = Iiri27N ..- 

WRI1E Li,212) FKI), ID, 

212 

FORMAT <E15‘: S, 2H ” , 10)- 

601 

IF < B J0G-1. 5)660, 660, 651 

IF L PLINC H• 0. 5 ) 650, t 50, 660 

660 

216 

IF L N- Nft.’C) 216, 216, 21 5 

IF' N - N ft! ) 1 j;, 21 '< j | -, 

r‘* 

NR-NR-t 1 


NRC-= 1 .~ - ~—- 

KONT =1 

IL N F - 4 .~ - —"— 

WRITE <1,201) NR 

215 

220 

CO 0 213 ~ - 

IFO ONT-ILNE) 214, 220. 21-4 

ilne-iLne+4 


WRITE <1, 221 > 


221 
21 4 


FORMfiT <77T 
RONT- KON T-tl 
NRC-NRC-tl “ 
CONTINUE 
11 = I1V1 


~ 2 a 0 

CONTINUF 

— „ ^ 

_GO TO < 14 5, 161) J OG 

‘ 65 8 

WRITE L 2, 6 6 5 ) ft < I ) 


_ GO TO 660 

620 

FORMfiT </-IS PUNCH 

625 

FORMAT < VES = 0N0 = 

630 

FORMAT < F2. 0) 

.63 5 

FORMAT < "PLACE DATA 

636 

FORMAT < TOP OF FOI 

__ 66 5 

FORMAT <E15 8) 


END " .~ 

# 

—‘-- -- . ___ 


PfiPER TftPE OUTPUT 


TAPE IN RE ftD £ ft 
•N' HIT CONTINUE- 


DES 

IO) 


I RED?'' ) 


































































